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ABSTRACT 

We present a panoramic study of the Fornax dwarf spheroidal galaxy, using data ob¬ 
tained as part of the VLT Survey Telescope (VST) ATLAS Survey. The data presented 
here - a subset of the full survey - uniformly cover a region of 25 square degrees centred 
on the galaxy, in r and i-bands. This large area coverage reveals two key differences 
to previous studies of Fornax. First, data extending beyond the nominal tidal radius of 
the dwarf highlight the presence of a second distinct red giant branch population. This 
bluer red giant branch appears to be coeval with the horizontal branch population. 
Second, a shell structure located approximately 1.4 degrees from the centre of Fornax 
is shown to be a mis-identihed background overdensity of galaxies. This last result 
casts further doubt on the hypothesis that Fornax underwent a gas-rich merger in its 
relatively recent past. 

Key words: galaxies: dwarf - galaxies: stellar content - Local Group - galaxies: 
individual (Fornax dSph) 


1 INTRODUCTION 

According to the standard ACDM cosmology, the Universe 
as we see it today is built via a scale-free process in which 
matter collapses under gravity and iteratively forms larger 
and larger structures. Dark matter dominates this hierar¬ 
chical structure formation process. Dwarf galaxies, the most 
dark matter dominated objects known, are therefore cru¬ 
cial laboratories for studying galaxy formation and evolution 
processes. Considerable work has been invested in the study 
of dwarf galaxies in the Local Group, summarised recently 
in 


Tolstoy, Hill V Tosi (2009) and McConnachie (2012). 


Local Group dwarf galaxies offer us the opportunity to 
examine resolved stellar populations in considerable detail. 
Such analyses have revealed richly complex star formation 
histories (e.g. Bono et al.|2010 Battaglia et al.|2011| de Boer 
et al.||2012[ to name just three recent examples), and allow 
us to search for evidence for tidal interactions and mergers. 
This last point is particularly crucial, since tidal disturbance 
on shallow potentials can directly impact on the interpre¬ 
tation of dwarf galaxy properties. Wide-held imaging from 
instruments such as the Dark Energy Camera (DECam, e.g. 
McMonigal et al.|2014 Roderick et al.|20T5 ) and survey tele¬ 


scopes such as the VLT Survey Telescope (VST) are making 
it possible to search for low surface brightness features in the 
outskirts of nearby dwarf galaxies. 

The Eornax dwarf spheroidal galaxy is one of the more 
intensively studied satellites of the Milky Way. Eirst discov¬ 


ered by Shapley (1938), it is relatively distant at 147±12kpc 


(all parameters in this paragraph from the review of |Mc-| 
Connachie 2012). Nevertheless, it is second only to the 


Sagittarius dwarf spheroidal in terms of luminosity {My — 
— 13.4 ± 0.3) and dynamical mass (Mdyn = 5.6 x lO^M© 
within the half light radius) amongst the MW dwarf popu¬ 
lation. Eor detailed information on the study of stellar pop¬ 
ulations in Eornax, we refer the reader to [Battaglia et ah] 
(2006), de Boer et al. (2012), and del Pino et al.| ( |2013 ). 

Eornax is of particular interest here because of a detec¬ 
tion of a possible large outer shell, located approximately 
1.4° from the centre of the dwarf (Coleman et al. 2005). 


The association of this shell with Eornax was circumstan¬ 
tially supported by an earlier detection of a smaller, inner 
feature with a shell-like morphology (hereafter referred to 
as clump 1, or Cl; Coleman et al. 2004). Shells are found 


E-mail: nbate@ast.cam.ac.uk (NFB) 


in large galaxies due to events such as mergers, and since 
hierarchical structure formation is scale free, we should ex¬ 
pect to see them in dwarf galaxies too. So far, Eornax is 
the only dwarf galaxy reported to display unambiguous ev- 
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Figure 1. Three-colour image of 16 square degrees centred on the Fornax dwarf spheroidal, built from VST ATLAS r and Tband 
imaging. The blue cross marks the centre of Fornax, and the blue ellipse marks the tidal radius ( [Battaglia et al.']|20Q6| . The blue box 
to the North-West marks the location of the outer shell identified in [Coleman et ar] ( |2005| >. This image was generated using 23 separate 
pointings from the VST ATLAS Survey. 


idence of a shell, although kinematic measurements of the 


Andromeda II dwarf galaxy (Amorisco, Evans <V van de Ven 


2014) suggest that it too may have undergone a relatively 


major merger. 

The initial interpretation of these two shells in Fornax 
was that the galaxy had undergone a gas-rich merger with 


another dwarf approximately 2 Gyr ago (Coleman et al. 


2004 Coleman et al. 2005). However, this hypothesis has 


proven difficult to reconcile with the dwarf galaxy’s star for¬ 
mation history, which shows a clear peak at intermediate 
ages; either 4 Gyr ago ( de Boer et al.[[2012 ), or 8 Gyr ago 
( del Pino et aL[|2013 ). 


Furthermore, Olszewski et al. (2006) demonstrated that 


the metallicity of the stars in Cl closely matched that of 
nearby field stars in Fornax. This result was confirmed by 


de Boer et al. (2013), who also detected a second overdensity 


in blue plume stars outside the centre of Fornax. This feature 
is composed of very young stars (~ 100 Myr old) with high 
metallicity, so is unlikely to have been formed from accreted 
pristine gas. 

A plausible explanation for these features, put forward 


Salvadori, Ferrara & Schneider (2008), is that these clumps 


were formed from enriched gas blown out of Fornax and 
subsequently re-accreted. Such gas would fall back into the 
dwarf over a timescale of ~ 250 Myr, sinking deeper as it is 
re-accreted. This is consistent with the radial stellar popu¬ 


lation age gradient observed by Battaglia et al. (2006), sug¬ 


gesting the gradual removal of gas from Fornax, however the 
timescales are difficult to reconcile with its complex star for¬ 
mation history. Disrupted clusters can also produce shells or 
clumps ( Penarrubia, Walker &: Gilmore|2009 ), however none 
of the globular clusters in Fornax show any signs of disrup¬ 


tion ( Mackey k. Gilmore||2003 Rodgers Roberts[[1994 ). 


Tidal interactions with the Milky Way are not expected 
to have played a significant role in the formation of any of 
these stellar substructures. After Sagittarius, Fornax is the 
Milky Way’s most massive dwarf spheroidal satellite. It is 
currently thought to be at or near perigalacticon, a large dis¬ 
tance 150 kpc, McGonnachie|2012 ) from the centre of the 
Galaxy, on a low-eccentricity polar orbit (e.g. Dinescu et al. 
2004 Piatek et aL|2007 ). There is, however, some conflict in 
proper motion measurements. Dinescu et al. (2004) results 


by Olszewski et al. (2006), Goleman <V de Jong (2008) and indicate that Fornax may have had an encounter with the 
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Magellanic Stream ~ 190 Myr ago, which could explain very 
recent bursts of star formation, and is supported by small 
scale structure found along the proposed orbit. However, 
these results are not consistent with other proper motion 
measurements (e.g. Piatek et al.||2002 ). 


Fornax’s strong population gradient (e.g. Battaglia 
et al.|200^ Held et al.|2010 Cesetti et al.|20lT ), with young 


stars centrally concentrated in the core, is not uncommon for 
dwarf spheroidals ( Harbeck et al.|20M ). Such gradients have 
been detected in the Local Group dwarf spheroidal popula¬ 
tion via imaging (e.g. [Harbeck et al.|2001 ) and spectroscopy 


le.g. 


Tolstoy et al.|2QQ4 ). Similarly, a clear metallicity gradi¬ 


ent is seen in Fornax (e.g. Walker, Mateo 01szewski|2009 


|de Boer et al.|20T2 ). Cesetti et al. (2011) detected bimodal¬ 
ity in the subgiant branch, and Battaglia et al. (2006) found 
a hint of splitting near the tip of the RGB. On the other 


hand, del Pino et al. (2013) and de Boer et al. (2012) find 


no clear evidence of splitting. 

In addition, Fornax is known to have distinct kinematic 
populations (e.g. Tolstoy et al.||2004 Battaglia et al.|[2006 


Walker Peharrubia||2011 Battaglia et al.||2(3TT ), with the 

lower velocity dispersion population being characterized by 
a lower metallicity. The major axes of these two populations 
are also offset by 30-40° ( Stetson et al.|1998 Battaglia et al. 


2006). This all points to a very complicated star formation 


history. 

The photometric studies discussed above have tended 


to focus on strategically placed sub-fields (e.g. de Boer 
et al.|20l2 Olszewski et al.|2006 ), or on coverage extending 


roughly to the Fornax tidal radius (e.g. [Battaglia et al.|2006 


Coleman de Jong|2008 ). In this paper, we use data taken 


on the 2.6-metre VLT Survey Telescope (VST) as part of the 
VST ATLAS Survey to explore a 25 square degree field cen¬ 
tred on Fornax. This area extends well beyond the nominal 
tidal radius, and beyond the ‘outer shell’ fields in jColemanj 


et al. (2005) and Coleman & de Jong (2008). It offers the 


perfect opportunity to cleanly explore the spatial distribu¬ 
tion of various Fornax stellar populations. 

In Section we discuss the VST data used throughout 
this paper, including the first re-analysis of the region sur¬ 
rounding the purported outer shell in Coleman et al. (2005), 


and discussions of colour magnitude diagrams and surface 
density profiles. In Sectionj^we use a matched filtering tech¬ 
nique to explore the spatial distribution of four populations: 
two distinct red giant branches populations revealed in our 
data, as well as the horizontal branch and young main se¬ 
quence. We use these blue plume stars to explore the inner 
portions of Fornax for overdensities in Sectionand discuss 
the implications of these measurements in Sectionj^ We con¬ 
clude in Section Throughout this paper, colour maps are 
generated using the ‘cubehelix’ scheme ( |Green|2011| . 


2 DATA 


The data presented in this paper were obtained as part of 
the ATLAS survey, one of three ESC public surveys cur¬ 
rently being carried out using the 2.6-m VLT Survey Tele¬ 
scope (VST). Details of the VST ATLAS survey and Data 


Release 1 (DRl) can be found in Shanks et al. (2013), Ko- 


posov et al. (2014) and Shanks et al. (2015). The survey 


aims to cover an area of several thousand square degrees in 
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Figure 2. Panel (a): a three colour image of the whole outer shell 
region identified in [Coleman et aL] ( |200^ . Panel (b): a zoom in 
on an example portion of the outer shell region, displaying an 
overdensity of galaxies. Panel (c): a field centred on the Fornax 
dwarf spheroidal, for comparison. The location of the (b) and (c) 
fields are shown in Figure [^ 
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(c) Central Fornax field 
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4 N. F. Bate et al. 


Table 1. Fornax VST ATLAS catalogue 


No. 

RA (ICRS) 

m 

Dec (ICRS) 

(O , „) 

9 

(mag) 

As 

(mag) 

g class 

r 

(mag) 

Ar 

(mag) 

r class 

i 

(mag) 

Ai 

(mag) 

i class 

527 

2 36 10.7473 

-31 59 56.741 

15.466 

0.001 

-1 

14.959 

0.001 

-1 

14.534 

0.001 

-1 

528 

2 33 20.4378 

-31 59 56.808 

22.050 

0.073 

1 

21.200 

0.062 

1 

20.719 

0.082 

1 

529 

2 34 22.3833 

-31 59 56.937 

21.597 

0.050 

-1 

20.164 

0.026 

-1 

19.108 

0.022 

-1 

530 

2 50 37.3408 

-31 59 57.123 

20.994 

0.023 

-1 

19.371 

0.013 

-1 

18.320 

0.011 

-1 

531 

2 32 19.9185 

-31 59 57.145 

19.285 

0.009 

-1 

18.003 

0.005 

-1 

17.291 

0.006 

-1 

532 

2 48 10.9424 

-31 59 57.187 

22.725 

0.102 

1 

21.612 

0.085 

1 

20.847 

0.098 

-3 

533 

2 36 34.4322 

-31 59 57.199 

22.845 

0.144 

-3 

21.311 

0.071 

-3 

20.477 

0.068 

-1 

534 

2 31 56.0799 

-31 59 57.236 

21.096 

0.033 

-1 

19.506 

0.016 

-1 

17.883 

0.009 

-1 

535 

2 28 18.9850 

-31 59 57.308 

0.000 

0.000 

0 

21.313 

0.068 

1 

20.315 

0.061 

1 

536 

2 35 19.1191 

-31 59 57.319 

20.389 

0.019 

1 

19.052 

0.011 

1 

18.286 

0.011 

1 


A sample of 10 entries in the final Fornax catalogue. The full catalogue is available as online supplementary material. 
Cambridge Astronomical Survey Unit pipeline classifications are described in Section]^ 


All objects Stars 



2 0 -2 2 0 -2 

^ (degrees) ^ (degrees) 


Figure 3. The spatial distribution of objects lying on the red giant branch and red clump locus in the Fornax field. Left panel: all 
catalogue detections. Right panel: all objects reliably classified by the CASU pipeline as stars. The ‘outer shell’ identified in |Coleman| 
|et al.| ( [2005| > is located ~ 1.4° north-west of the centre of Fornax in the left panel. It does not appear in the right panel. Red boxes mark 
the location of the black and white zoom fields provided in Figure 


the southern celestial hemisphere, in ugriz filters to depths 
comparable to the Sloan Digital Sky Survey (SDSS) in the 
north. The median seeing in DRl is < and the median 
limiting magnitudes in each band corresponding to 5cr de¬ 
tection limits are approximately 21.0, 23.1, 22.4, 21.4, 20.2. 


The raw data are processed by the Cambridge Astro¬ 
nomical Survey Unit (CASU; [lrwin &: Lewis|200T ). Objects 
are parameterised and classified morphologically, and cat¬ 
alogues are generated. The ATLAS DRl photometry used 
throughout this paper is in Vega magnitudes. These magni¬ 
tudes were corrected for Galactic extinction using [Schlegel,| 
Finkbeiner <V Davis (1998) dust maps in conjunction with 


SDSS filter extinction coefficients obtained from |Schlafiy <V| 
Finkbein^ ( |2011 ). 


Table contains a sample of 10 catalogue entries. The 


full catalogue is available as online supplementary mate¬ 
rial. For each catalogue object we provide ICRS coordinates, 
gri photometry and associated root-mean-square errors, and 
CASU pipeline classification in each of the three bands. The 
pipeline classifications are: 1 = galaxy; 0 = noise-like; — 1 
or — 2 = reliable star; — 3 = star/compact galaxy; — 6 = low 
average confidence value; —7 = contains bad pixels; —8 = 
greater than l" matching error (matching criteria search ra¬ 
dius is 2.5''); —9 = saturated. These supplementary data are 
supplied without correction for Galactic extinction. 


2.1 The Fornax outer shell 

In Figure^ we present a three-colour image constructed from 
the g, r and i filters of a 16 square degree region centred 
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Table 2. Fornax strutural parameters 


Parameter 


King 


King 

(corrected) 


Sersic 


Sersic 

(corrected) 


Battaglia et al. (20061 


O 2000 

2^39^52^ 

2h39"^51® 

2^39^52^ 

2^39^52^ 

2^39^52^ 

(^2000 

-34°30'39" 

-34°30'39" 

-34°30'38" 

-34°30'38" 

-34°30'49" 

Ellipticity 

0.31 + 0.01 

0.31 + 0.01 

0.31 + 0.01 

0.31 + 0.01 

0.30 + 0.01 

Position angle 

41.6° ±0.2° 

41.5° +0.2° 

41.5° +0.2° 

41.6° +0.2° 

46.8° + 1.6° 

Tc 

15.4' ±0.1' 

14.6' + 0.1' 

- 

- 

17.6' + 0.2' 

n 

69.5' ± 0.4' 

69.7' + 0.3' 

- 

- 

69.1' + 0.4' 

Fs 

- 

- 

15.2' + 0.1' 

14.5' + 0.1' 

17.3' + 0.2' 

Sersic index m 

- 

- 

0.760 + 0.004 

0.779 + 0.005 

0.71 + 0.01 

Contamination surface density 

0.501 + 0.001 

0.501 + 0.001 

0.499 + 0.001 

0.499 + 0.001 

0.78 + 0.02 


(stars arcmin ^) 


Structural parameters obtained from fits to the photometric catalogue above the VST ATLAS 5cr detection limits. In the 
corrected case, a modification is made to the models to account for crowding in the inner regions of Fornax (see Section [2.2| 
for details). 


on the Fornax dwarf spheroidal. The coordinates (^, 77 ) rep¬ 
resent a tangent plane projection centred on Fornax. This 
image was stitched together from 23 separate pointings in 
the VST/ATLAS Survey. Of particular interest in this wide- 
held image is the ‘outer shell’ first noted in |Coleman et ah] 
( 20Q5| ; its position is marked with a blue box to the North- 
West of the Fornax dwarf. The dwarf spheroidal itself is 
clearly visible in this image, and there is perhaps some hint 
of an overdensity in the blue boxed region. In Figure we 
zoom in on this region. 


The CASU pipeline uses a curve of growth analysis to 
classify objects as noise detections, galaxies, and probable 
stars. To reflect the method used by Coleman et al. (2005) in 
their original detection of the Fornax outer shell, in Figure 
1 ^ (left panel) we plot all catalogue detections that lie on the 
Fornax red giant branch (RGB) and red clump (RC) locus 
in the g and r filters (see Figure]^ for the relevant colours 


and magnitudes). The feature identified in Coleman et al. 
( 2005| as a shell associated with Fornax is clearly visible as 
an extended feature approximately 1.4° to the North-West 
of Fornax. 


In the right panel of Figurewe instead plot all objects 
lying on the same RGB+RC locus that were reliably clas¬ 
sified by the CASU pipeline as stars (aperture photometry 
classifications of —1 or —2 in g, r and z, which corresponds 
to point sources up to 2cr from the stellar locus). It is im¬ 
mediately obvious that the ‘outer shell’ has disappeared. 


We therefore conclude that the Fornax outer shell does 
not exist; it is in fact a background over density of galaxies. 
This casts doubt on the possibility of a major merger in the 
past of the Fornax dwarf spheroidal. The implications of this 
finding will be discussed further in Section 

The result presented in Figure does not depend on 
the magnitude cuts chosen for the RGB+RC. The cut cho¬ 
sen in that figure reaches below the red clump, to g = 22.3 
magnitudes. If we instead cut above the red clump, select¬ 
ing only the red giant branch stars {g < 21.1 magnitudes), 
the features in Figure are preserved. As a final check, we 
visually inspected the VST ATLAS images taken in the best 
seeing. An example image is provided in Figure [2(1^ where 
an overdensity of galaxies is clearly visible; for comparison 
a field centred on Fornax is provided in Figure |2(c) This 


inspection confirms the CASU pipeline identifications in the 
vicinity of the ‘outer shell’. 

For the remainder of our analysis, we will use all objects 
in the final catalogue reliably classified as stars in all three 
bands. We note again that these are aperture photometry 
classifications of —1 or —2 in g, r and z, which corresponds 
to point sources up to 2 cr from the stellar locus. 


2.2 Surface density profiles 

We fit a surface density profile to the data using an affine- 
invariant ensemble sampler for Markov chain Monte Carlo 


(MCMC) of the form proposed by Goodman <V Weare 
( |201Q| ). This allowed us to easily sample the entire parame¬ 
ter space, while leaving free the parameters for the centre of 
Fornax (right ascension and declination), the position angle, 
the ellipticity, the contamination stellar density, and three 
parameters for the radial density profile - a total of 8 param¬ 
eters. To get a good fit, we ran the MCMC with an ensemble 
of 50 walkers, each of which did approximately 1000 steps. 

In line with Battaglia et al. (2006), we fit the Fornax 
stellar densities using both a King profile ( King|19^ ) and a 
Sersic profile ( Sersic|19^ ). For a detailed explanation of the 
reasons behind these choices, we refer the reader to |Battaglia| 
et al. (2006). We define the ellipticity e to be e = 1 — 5/a, 
where h and a are the semi-minor and semi-major axes of 
the galaxy respectively. The position angle is taken from 
North towards East (counter-clockwise, in the figures pre¬ 
sented throughout this paper). 

Following the method outlined in the appendix of 
Richardson et al. ( 2011| , we calculated likelihoods in a way 
which does not require binning, thus using the full set of 
information available from the data within the 25 square 
degree region around Fornax. The only cuts performed were 
to select for entries classified as stellar (—1 or — 2 ) in all 
three filters, and to exclude entries which fell below the 5 cr 
detection levels of the VST survey in any of the three filters 
(g = 23.1, r = 22.4, z = 21.4). The fit was also performed 
using only the stellar classification cuts, and produced con¬ 
sistent values for all of the parameters with the exception 
of the contamination stellar density, which was 2 per cent 
higher in this case. 

To account for crowding near the core, we applied the 
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Figure 4. Best King and Sersic profile fits, compared with the data binned in elliptical annuli using the best fit position angle and 
ellipticity. The fits themselves were performed to the full unbinned dataset above the VST ATLAS 5a detection limits, allowing all 
parameters to vary freely, as described in Section [2^ Fits are provided for models both with and without a correction for crowding in 
the centre of Fornax, although the correction is minor. Parameters for the fits are provided in Table 


correction detailed in the appendix of [Irwin fc Trimble] 
(1984) to the models and repeated all the fits. Again we 


found that most parameters remained consistent for all fits, 
with a 2 per cent increase in the contamination in the ab¬ 
sence of any magnitude cuts. For the King profile, applying 
the crowding correction results in a slightly smaller core ra¬ 
dius Fc- In the Sersic case, the crowding correction leads to a 
correspondingly smaller scale radius and a 2a increase in 
the Sersic index m. This is unsurprising, since the crowding 
correction effectively increases the number of stars in the 
Fornax core. 

We present relevant parameters from the King and Ser¬ 
sic hts in Table where the conhdence intervals represent 
the error in the htting procedure. The results for the ap¬ 


propriate models from Battaglia et al. (2006) are also pro¬ 


vided for comparison. Our results are broadly consistent 
with that work, whose coverage did not extend uniformly 
beyond the Fornax tidal radius. We note in particular that 
there is a ~ 5° offset between our measurements of posi¬ 
tion angle. Our results are consistent with the earlier work 


of Irwin & Hatzidimitriou ([1995 ), rather than the hgures 


quoted in Mateo (1998) and [Battaglia et al. (2006). We also 


hnd slightly smaller King core radius and Sersic scale radius 
than Battaglia et al. ( 2006| . 

For visual comparison. Figureshows the best ht mod¬ 
els for each prohle compared to the data binned in ellipti¬ 
cal annuli. These hts were obtained by applying magnitude 
cuts at the VST ATLAS 5a detection limits, and are shown 
both with and without the crowding correction. As noted in 
Battaglia et al. (2006), it is clear that the Sersic prohle pro¬ 


vides a much better ht around the tidal radius, whereas the 
King prohle is a slightly better ht at small radii, although 
the data show a hatter core region than either prohle can 
accommodate. 


2.3 CMDs 

In Figures we provide colour magnitude diagrams 

(CMDs) in three hlter combinations: g and r, r and z, g and 
i respectively. Each hgure shows three CMDs: for the central 
regions of Fornax (elliptical radius Veii < 0.20°, left panels), 
an annulus extending to the tidal radius (0.40° < Veii < 
1.15°, middle panels), and a contamination held consisting 
of ah data with ^ < —2.00° (right panels). The contamina¬ 
tion held was chosen to best match the depth of the central 
Fornax helds. We note that while the spatial area covered 
by the CMDs in the second and third panel is roughly the 
same (2.56 and 2.50 square degrees), the hrst panel, centred 
on Fornax, corresponds to a much smaller area (0.09 square 
degrees). 

The division into a central region and an outer elliptical 
annulus was chosen to highlight structure in the red giant 
branch (RGB). In Figuresand(both p-band), the outer 
CMD shows a clear bifurcation in the RGB. There were hints 
of this feature in previous work (e.g. [Battaglia et al.[[2006 


Coleman <V de Jong[[2008 ), however this is the hrst time it 


has been revealed with such clarity. Indeed, previous work 
covering smaller spatial regions 


(.e.g. 


Stetson et al. 1998 


de Boer et al.|20'T2 ) has favoured continuous star formation 


at early ages, although de Boer et al. (2012) note that this 
could be a consequence of the age resolution in their star 
formation histories. 


Following Battaglia et al. (2006), we will refer to the two 


RGBs as the blue RGB (B-RGB) and red RGB (R-RGB). 
The former is less dense and distributed over a wider area, 
as expected for an older/more metal poor population. Con¬ 
versely, although the R-RGB population is visible at ah el¬ 
liptical radii, it is much more centrally concentrated. Young 
blue plume stars are found only in the innermost regions 
of Fornax, whereas the horizontal branch is virtually absent 
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< 0.20° 0.40° < < 1.15° f < -2.00° 



0 1 2 0 1 2 0 1 2 
g-r g-r g-r 


Figure 5. Colour magnitude diagrams in g and r for two regions around Fornax, and a contamination region away from the dwarf. Left 
panel: all stars within an elliptical radius of 0.20 degrees centred on Fornax. Middle panel: all stars within an annulus with inner elliptical 
radius 0.40 degrees, and an outer elliptical radius equal to the Fornax tidal radius {rt = 69.1^ [Battaglia et al.||2006| >. Right panel: all 
stars in a rectangle with —2.50 < ^ < —2.00 and —2.50 < 77 < 2.50. 
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Figure 6. Colour magnitude diagrams in r and z, as per Figure 
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Figure 7. Colour magnitude diagrams in g and z, as per Figure 
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Figure 8. Colour magnitude diagrams for two regions around Fornax. Left panels: all stars within an elliptical radius of 0.20 degrees 
centred on Fornax. Right panels: all stars within an annulus with inner elliptical radius 0.40 degrees, and an outer elliptical radius equal 
to the Fornax tidal radius {rt = 69.1^ [Battaglia et al.|2006| >. A variety of features are marked. HB: horizontal branch, BP: blue plume (a 
young main sequence, not to be confused with the blue stragglers at the end of the older main sequence turn off), BL: blue loop Helium 
burning stars, RC: the red clump, AGBb: asymptotic giant branch bump, B-RGB: bluer red giant branch, R-RGB: redder red giant 
branch (the red giant branch continues below the red clump). In the outer annulus, the metal richer R-RGB population is diminished 
and a second metal-poorer B-RGB is visible. The bottom row contains Hess diagrams, constructed with 0.01 magnitude pixels and 
square-root scaling, smoothed using a Gaussian kernel with characteristic radius 0.04 magnitudes. No contamination has been subtracted 
from these Hess diagrams. 


in these same regions. These features also point to younger 
stars being preferentially found in the centre of Fornax, as 
expected. 

In Figurewe again provide g — r CMDs for the central 
region and tidal radius annulus (top panels), along with Hess 
diagrams of the same (bottom panels). The square root scal¬ 
ing of the colour gradient on the Hess diagrams was again 
chosen to highlight the bifurcation in the RGB, which is 
most pronounced in this combination of hlters. 

In the top panels of Figure we highlight a variety of 
CMD features, including a blue main sequence of young stars 
(referred to hereafter as blue plume stars) and an extended 


blue horizontal branch indicative of an ancient stellar popu¬ 
lation. These features will be explored in more detail in the 
following section, wherein we use the Hess diagrams in the 
bottom panels of Figure to construct matched hlters. 


3 SPATIAL DISTRIBUTION 

To explore the spatial distribution of the B-RGB, R-RGB, 
horizontal branch, and blue plume Fornax populations in 
a statistically robust way, we perform a Poisson-based 
matched hltering process. This process was hrst described 
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McMonigal et al. (2014); the salient details are repeated 


here. 


3.1 Matched filtering technique 


We wish to determine the spatial extent of various popula¬ 
tions visible in the Fornax CMDs in Section [2.3| To do so, 
we apply a variation on the matched filtering technique of 


Kepner et al. (1999) and Rockosi et al. (2002). 


Matched filtering in this context is a technique where 
observed data for an object is used to construct a filter that 
describes the shape of the population of interest in colour- 
magnitude space. This filter, in combination with an appro¬ 
priate contamination, is applied to low-signal regions to pick 
out the source population from the noise in a statistically 
meaningful way. The standard matched filtering technique 
as described in Rockosi et al. ( 2002 ) and Odenkirchen et al. 


(2003) assumes Gaussian statistics, however our CMDs are 


governed by Poissonian processes. We therefore modify the 
matched filtering technique as follows. 

We describe the number of stars Xc,m in a given (^, 77 ) 
spatial pixel, with a given colour and magnitude (c, m) with 
the following model: 


Xc,m,^,r] — T ^c,7i 


( 1 ) 


iT,m is the foreground filter, in the form of a probability 
density function in colour-magnitude space. The contamina¬ 
tion is described by and is the number of source 

stars in the spatial pixel It is this last, that we 

wish to determine. 

The algorithm we use to determine is: 

(i) For a given spatial pixel (^, 77 ), construct an observed 
Hess diagram. 

(ii) Set a^^r] = 0. 

(hi) Calculate Xc,m,^,r] according to Equation 

(iv) Determine the probability for the observed number 
of stars given a^^r] in each colour-magnitude pixel. 

(v) Calculate the probability for a^^r] over the entire 
colour-magnitude space. 

(vi) Increment by 1, and return to step (ii) until a 
maximum is reached (~ 1000). 

We calculate the probability for the observed number 
of stars nc,m,i,r] given a^^r] using a Poisson distribution: 

P(nc,m,^,77 |o^^,T7) = exp( Xc^m,^,ri) • (2) 

The probability for in a given spatial pixel is the 
product of the individual probabilities at each pixel in colour 
magnitude space: 

p(jlc,rri,i,r)\0L^^ri) ■ (3) 

We take the maximum likelihood value of around which 
the probability distributions are sharply peaked, to be the 
number of source stars in a given spatial pixel. 


3.2 Filters 

Filters were generated from {g — r, g) Hess diagrams, as 
displayed in Figure and the top panels of Figure This 



g-r g-r g-r 

Inner-outer (R-RGB) Outer-inner (B-RGB) 


03 



^ 0.5 1 1.5 0.5 1 1.5 


g-r g-r 

Figure 9. Hess diagrams illustrating the process used to de¬ 
termine appropriate regions to use for red RGB (R-RGB) and 
blue RGB (B-RGB) filters. Top panels show Hess diagrams con¬ 
structed from the innermost region of Fornax {top left), an outer 
elliptical annulus around Fornax {top middle), and a contamina¬ 
tion region north of Fornax {top right). Filters were generated 
using these Hess diagrams. Bottom panels show subtractions of 
appropriately scaled Hess diagrams to aid in the selection of the 
filter regions, marked in white. We note that the bottom panels 
are illustrative only; the filters themselves were constructed from 
the original Hess diagrams in the top panels. 


combination of filters was found to best discriminate be¬ 
tween the two populations in the split RGB. We repeated 
the matched filtering process for the {g — i,i) case, and found 
comparable results. 

Hess diagrams were constructed for three spatial regions 
in the data: the innermost region {reii < 0 . 20 °), an annulus 
out to the Fornax tidal radius (0.40° < reii < 1.15°), and a 
contamination field away from the Fornax population (^ < 
—2.00°). These Hess diagrams consisted of 0.025 x 0.025 
mag pixels in both colour and magnitude, smoothed using 
a Gaussian with 2 pixel dispersion. This was comparable to 
the observational errors. 

The filters were built by truncating the appropriate 
Hess diagram at appropriate GMD boxes. Figure illus¬ 
trates how these boxes were chosen for the B-RGB and R- 
RGB populations. By using scaled subtractions of the inner 
ellipse and outer annulus Hess diagrams (bottom panels), 
we were able to select regions where one population or the 
other dominated. The final R-RGB and blue plume filters 
were built from the (raw) inner Hess diagram, and the final 
B-RGB horizontal branch filters were built from the (raw) 
outer annulus Hess diagram. The regions selected for the 
horizontal branch and blue plume populations are illustrated 
in Figure 
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Figure 10. Top left panel: matched filtered stellar density contours for the blue red giant branch (B-RGB) population. Top right panel: 
the red red giant branch (R-RGB). Bottom left panel: the horizontal branch (HB). Bottom right panel: the young main sequence (blue 
plume, BP). Pixels are 2' x 2', smoothed using a Gaussian kernel with 4' dispersion. Gontours are displayed at 2, 5, 10, 25, 50, 100, 200, 
300 times the (smoothed) RMS value calculated in the contamination region (see text). The tidal ellipse from [Battaglia et al.| ( [2006| > 
is marked with a dashed red line, and from this paper with a solid blue line. The centroid (red cross) is marked at the parameters 
determined in this paper. The known Fornax globular clusters are marked with filled blue circles. 


Foreground filters are probability density functions, and 
so once a region has been selected from a Hess diagram, 
they are normalised independently so the area underneath 
the filter is equal to 1. Conversely, the contamination filter 
represents a number of stars that is expected in each colour 
and magnitude pixel, and so it is simply normalised to unit 
area. 


We attempted to push as far down the red giant branch 
as possible, to maximise our signal while hopefully avoid¬ 
ing contamination from the red clump. As a test, we re¬ 
peated our matched filtering analysis using filters truncated 
at g = 20.5 (i.e. removing the vertical section of the selection 
regions). Although the signal was slightly diminished in the 


resulting matched filtered maps, the structure of the maps 
was essentially unchanged. 


3.3 Matched filtered maps 

Matched filtered maps for four populations are provided in 
Figure [To) the B-RGB (older/more metal-poor stars; top left 
panel), R-RGB (younger/more metal-rich stars, top right 
panel), horizontal branch (indicative of an ancient stellar 
population, bottom left panel), and the blue plume (a young 
main sequence, bottom right panel). We note that the pan¬ 
els are normalised separately, based on their respective root- 
mean-square contamination in the ^ < —2.00° contamina¬ 
tion region. This level is 0.072 stars pixel”^ for the B-RGB 
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Figure 11. Stellar number count profiles for the blue red giant branch (B-RGB, left panel) and red red giant branch (R-RGB, right 
panel) Fornax populations, as defined in Section [3.2| These profiles were obtained from matched filtered maps, in a 20' slice along the 
major axis. The dashed lines mark the RMS stellar density in the contamination region of each map. These RMS values were used to 
determine the contours displayed in Figure 



Elliptical semi—major axis radius (arcmin) 

Figure 12. The cumulative fraction of stars as a function of el¬ 
liptical semi-major axis radius for four populations: the blue red 
giant branch (B-RGB, black squares), the red red giant branch 
(R-RGB, red circles), the blue plume (BP, purple crosses), and 
the horizontal branch (HB, blue triangles). Number counts were 
obtained from matched filtered maps, accounting for the expected 
RMS contamination in each ellipse. Each population is normalised 
independently. We note that the B-RGB and HB follow approxi¬ 
mately the same profile, and so are likely to be coeval. 


(top left), 0.084 stars pixeP^ for the R-RGB (top right), 
0.034 stars pixel”^ for the horizontal branch (bottom left), 
and 0.026 stars pixel”^ for the blue plume (bottom right). 
In each panel we also plot the tidal ellipse and centroid de¬ 
termined in this paper, and the five known Fornax globular 
clusters. 

It is clear from these maps the the B-RGB population 
is significantly more diffuse than the R-RGB population. 
The R-RGB population is more centrally located, and dis¬ 
plays somewhat less regular density contours, particularly 
towards the outer regions. This suggests somewhat patchy 
star formation. Neither population shows any clear evidence 


of extra-tidal features, suggesting that Fornax is not cur¬ 
rently undergoing tidal disruption. This is in line with proper 
motion measurements placing Fornax at a relatively distant 
perigalacticon of ~ 150kpc (e.g. Dinescu et al.|2004 ). 

The horizontal branch, while consisting of fewer stars 
than the B-RGB, appears to follow essentially the same spa¬ 
tial distribution as that population. This suggests that the 
two populations are coeval. The blue plume stars are once 
again quite centrally located, and their density contours sig¬ 
nificantly more circular than any of the other populations 
examined here. 

In Figure m we provide stellar radial profiles for the 
red giant branch populations, obtained from the matched fil¬ 
tered maps. These profiles were determined along a 20^-wide 
slice through the matched filtered maps, along the major 
axis. The left panel shows the B-RGB population, and the 
right panel the R-RGB population. The latter is denser and 
more centrally located, and appears slightly asymmetric, ex¬ 
tending out to the tidal radius on towards the South-West, 
but truncating earlier on the North-East. 

We also map the cumulative number of stars in each of 
the four populations, normalised to their peak number, as a 
function of elliptical semi-major axis radius in 2' increments. 
These prohles are shown in Figure They highlight that 
the blue plume population is the most compact (although we 
note again that it does not have the same elliptical shape as 
the other three populations), followed by the metal-richer R- 
RGB. The metal-poorer B-RGB and the horizontal branch 
follow essentially the same prohle as a function of elliptical 
radius, suggesting that these two populations are coeval. 


4 INNER OVERDENSITIES 


In Section 12.11 we discussed the misclassihcation of a back¬ 
ground overdensity of galaxies as a shell of stars, presumed 
to have been associated with a signihcant merger in For¬ 


nax’s past. Both Goleman et al. (2004) and de Boer et al. 
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Figure 13. The spatial distribution of Fornax blue plume stars, selected from the g and r bands. Top panels show individual catalogue 
objects. Bottom panels, the stars have been binned into 1.5' x 1.5' pixels, and a linear colour gradient chosen to highlight overdensities. 
Left panels: the entire blue plume population (excluding the horizontal branch). Middle panels: blue plume fainter than the horizontal 
branch. Right panels: blue plume brighter than the horizontal branch. In each panel, an ellipse is drawn with Fornax’s ellipticity and 
position angle (this paper), at half its measured tidal radius. Filled blue circles in the top right panel mark the locations of known 
globular clusters. Four overdensities are highlighted with 3' diameter circles, and labelled in the bottom right panel. Overdensity 1 was 
first detected in [Coleman et aT] (|20Q4| >, overdensity 2 in jde Boer et al.j ( [^Q13| >. Overdensities 3 and 4 are new. CMDs for these over dense 
regions can be found in Figure |14| 


(2013) have reported the detection of stellar overdensities 
much closer to the centre of Fornax. In the case of ide Bo"^ 
let al.| ( |2013l l, a new overdensity was identified in blue plume 
(BP) stars. 


Guided by this discovery, we examined the spatial dis¬ 
tribution of the blue plume population. CMD selection boxes 
were built using the central square degree of our g and r data 
(see Figure 14 top left panel), where the BP population is 
most easily discriminated. We split the BP into three popu¬ 
lations: the entire BP (excluding the horizontal branch), the 
BP below the horizontal branch (“faint blue plume”), and 
the BP above the horizontal branch (“bright blue plume”). 
We were motivated to look at the brighter and fainter BP 
populations separately by the unusually bright blue plume 
in the de Boer et al. (2013) overdensity when compared with 
the rest of the Fornax population. 


In Figure we present maps of the blue plume stars in 
these three CMD regions: the full blue plume (left panels), 
the faint blue plume (middle panels), and the bright blue 
plume (right panels). In the top panels, the individual stars 


are plotted, whereas in the bottom panels they are binned 
in 1.5^ X 1.5^ pixels, with a linear colour gradient chosen to 
highlight regions of higher density. In each panel, an ellipse 
is plotted with the ellipticity and position angle of Fornax, 
at an elliptical radius of half the Fornax tidal radius. 


Both the Coleman et ar] ( 20Q4| ) and de Boer et al. (2013) 
overdensities (circles I and 2 respectively) are clearly visible 
in the full BP population. The Coleman et al. ( 2004| region 
particularly stands out in the fainter BP, and the jde Boerj 


et al. (2013) region is a marked overdensity in the brighter 


BP population. This is in agreement with jde Boer et al.j 
(2013), who showed that while their new overdensity con¬ 
sisted of unusually young stars, the Coleman et al. (2004) 
overdensity was much more similar to the standard Fornax 
population. 


The Coleman et al. (2004) feature does appear to have 
a shell-like morphology, whereas the de Boer et al. (2013) 
feature is more radial, although still extended. Chip gaps 
certainly have an impact on structure we see in the VST 
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Figure 14. CMDs for overdensities in the Fornax field, following the labelling convention in Figure Top left: the full CMD for the 
central square degree around Fornax. The blue dashed box marks the blue plume (BP), and the solid red box the horizontal branch (HB). 
The HB div ides the blue plume in to two regions: the Brighter BP and the Fainter BP. Top mi ddle: a 3^ radius arou nd the overdensity 
identified in [Coleman et aL| ( |2004| ). Top right: a 3' radius around the overdensity identified in jde Boer et al!| ( |2013| >. Bottom left: a 3' 
radius around the overdensity in Figure [^identified in the Brighter BP population. Bottom middle: a 3' radius around the overdensity 
in Figure [T^ identified in the Fainter BP population. Bottom right: a 3' radius around the Fornax centroid identified in this paper. Not 
the similarities between overdensities 2 and 3, 1 and 4. 


ATLAS data, although the binned stellar density plots in 
Figure [Ts] do indicate a patchy distribution of young stars. 

We identify two new regions away from the very centre 
of Fornax that appear over-dense, one in the brighter BP 
population (circle 3) and one in the fainter BP population 
(circle 4). Neither is as significant as the earlier detections. 
The feature in circle 3, in the bright BP, is more diffuse than 
any of the other overdensities. The feature in circle 4 does 
appear to have a somewhat arc-like morphology, although 
the low number count of stars could well be misleading. 

In Figurewe display CMDs for the central square de¬ 
gree of Fornax, and for 3' circles centred on each of the over¬ 


densities in Figure [Ts] The bright blue plume in the jde Boerj 
et al. ( 2013| overdensity is striking (top right panel). Our 


new overdensity 3 (bottom left panel) is somewhat similar, 
suggesting that very recent star formation may be occur¬ 
ring in more than one isolated pocket outside the centre of 
Fornax. The least significant of the new overdensities that 
we’ve highlighted here, number 4 (bottom middle panel), is 
very similar to the feature hrst identihed in [Coleman et al.j 
([20041. 
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5 DISCUSSION 

The most prominent result of our analysis of wide-held VST 
ATLAS data of the Fornax dwarf spheroidal concerns the 


pair. This counter-rotating result hinges on the proper mo¬ 


tion measurements of Piatek et al. (2007). Clearly, more 


outer shell hrst reported in Coleman et al. (2005). With 


these new data, taken on a larger telescope in better ob¬ 
serving conditions, the outer ‘shell’ is revealed to be a back¬ 
ground overdensity of galaxies (see Figure]^, and not as¬ 
sociated with Fornax at all. This putative outer shell had 
been taken as evidence for a major interaction in Fornax’s 
past, and the gold-star example of dwarf-dwarf interactions 


in the Local Group (e.g. Deason, Wetzel V Garrison-Kimmel 


20141 along with Andromeda II[Amorisco, Evans V van de 
Ven|2014| |. 


A key feature of the Fornax merger hypothesis put for¬ 


ward by Coleman et al. (2004) and Coleman et al. (2005) 


was the presence of shell-like structures on opposite sides of 
the Fornax minor axis: a smaller inner shell on the southeast, 
and the much more prominent outer ‘shell’ to the northwest. 
These shells were also associated with ‘lobes’, consisting of 
slight overdensities of stars outside the tidal ellipse along 


spectroscopic data are needed, focussing in particular on 
the more metal poor red giants located in the outskirts of 
Fornax, where spectroscopic coverage tends to be sparse. 

Nevertheless, by demonstrating that the Fornax outer 
shell does not in fact exist, we cast doubt on the possibil¬ 
ity of a significant merging event in the history of Fornax. 
The causes of the young bursts of star formation and the 
(possible) counter-rotation in the metal poor and interme¬ 
diate metallicity red giants are still unclear, although the 
weight of evidence supports the conclusion that all of these 
stars were generated from Fornax gas, rather than material 
accreted from elsewhere. 

The wide-held VST ATLAS data that we present also 
clearly reveals a bifurcation in the Fornax red giant branch. 
There have been hints of this division in previous work (e.g. 
Battaglia et al.|[2006 Walker V Peharrubia||2Qll ), however 


it is clearly visible in the VST {g — r^g) CMDs. The clar¬ 
ity of this feature is primarily due to the uniform wide-held 


axis. Although the ESO/MPG 2.2m observations 

displays a radial age gradient (e.g. Stetson et al. [1998 [Sa- 

1 Coleman V de Jong| (|2008| specihcally cover the 

viane. Held & Bertelli[[2000[ [Battaglia et al.[[2006[ [de Boer 

1’, no information on that feature was reported. 

et al.[[2012[ [del Pino et al.|2013). 


by Yozin V Bekki (2012), in which they use N-body simu¬ 
lations to conclude that the hypothesis is plausible. This 
result hinged almost entirely on the creation of an outer 
shell, and further predicted the presence of tidal tails. We 
hnd no clear evidence for tidal features in our analysis, in 
line with previous observational and theoretical studies (see 
e.g. Peharrubia et al.|]2009 ). 

It should be noted that even despite the incorrect identi- 
hcation of the outer ‘shell’, reconciling a merger event with 
the Fornax star formation histories and the properties of 
the inner shells has proven problematic. [Olszewski et al^ 
( 20Q6| demonstrated that the stars in the inner shell identi¬ 
fied by [Coleman et aL (2004) were consistent with younger 
(~ 1.4Gyr) Fornax field stars, suggesting that they were the 
product of pre-enriched gas in Fornax rather than externally 
accreted. 


(2012) and del Pino et al. (2013) differ in their details. 


de Boer et al. (2012) prefer an intermediate age of 1-10 Gyr 


for the dominant Fornax population, with a dominant RGB 
population at ~ 4Gyr. They do note, however, that their 
age resolution may smooth out any old features in the star 
formation history. 


In contrast, del Pino et al. (2013) prefer a peak in star 


formation at approximately 8 Gyr ago in the central regions 
of Fornax, and 10 Gyr and older in outer regions (located 
roughly at the Fornax tidal radius). They attribute this dif¬ 
ference to shallower depth in the [de Boer et al.' ( [2012| pho¬ 
tometry. We note that our uniform, wide-held VST data only 
clearly reveal the second distinct peak in the RGB when the 
data are clearly extended to the tidal radius; repeating the 


analysis of del Pino et al. (2013) all the way out to that ra¬ 


dius may well heighten the signal in the older population. It 


This result is conhrmed in de Boer et al. (2013), where 


is perhaps signihcant to note that the del Pino et al. (2013) 


the authors conduct a detailed analysis of the CMDs and 
star formation histories of the inner Coleman et al. (2004) 


clump. They also hnd a new clump, which is even younger 
still (100-200 Myr, the youngest age bin considered in their 
analysis). By hltering only for bright blue plume stars, we 
reveal here a third clump of young star formation, ohset 
by ~ 7' from the centre of Fornax. This feature displays a 


CMD very similar to the de Boer et al. (2013) overdensity 
(see Figure p^. 

Although it is tempting to discard the merger hypoth¬ 
esis altogether, the situation is complicated somewhat by 
kinematic data. Using a technique developed by [Walker ^ 
Peharrubia (2011), Amorisco V Evans (2012) presented ev¬ 


idence that a sample of Fornax red giants consists of three 
distinct subpopulations, each with diherent kinematic prop- 


results, while not ruling out mergers as signihcant in the 
evolution of Fornax, are completely consistent with evolu¬ 
tion where mergers did not play a signihcant role. 


6 CONCLUSION 

We present photometric data from the VLT Survey Tele¬ 
scope ATLAS Survey of a 25 square degree region centred 
on the Fornax dwarf spheroidal. These uniform, wide-held 
data reach down to depths of 23.1 in 22.4 in r, and 21.4 
in i, and reveal a variety of interesting features. 

Most signihcantly, we demonstrate that the outer ‘shell’ 


erties. In particular, in two separate datasets (Walker, Mateo 
fc Olszewski[[20Q9[ and Battaglia et al.[[2006 ) they hnd that 


two red giant populations (one metal poor and one interme¬ 
diate metallicity) are rotating in opposite directions. They 
interpret this as possible evidence for the merger of a bound 


hrst reported in Coleman et al. (2005) is in fact not a stellar 
overdensity, but an overdensity of background galaxies. Most 
likely the misidentihcation in the original data was due to 
the observing conditions. No follow-up has been reported on 
this feature since its original appearance in the literature, 
although its presumed existence has been used to support 
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arguments that Fornax underwent a significant merger with 
another dwarf-sized object in its past. 

Multiple lines of evidence are suggestive of some sort of 
merger activity, including isolated regions of young star for¬ 
mation away from the core of Fornax ( Coleman et al.|2004 


de Boer et al. Pol'S this work), and the possible detection of 


counter-rotation in the metal poor and intermediate metal- 
licity populations ( Amorisco Evans||2012 ). Although our 
new results do not rule out the merger hypothesis entirely, 
they certainly remove a line of supporting evidence. 

We also detect for the first time a completely unam¬ 
biguous bifurcation in the red giant branch. This suggests 
that rather than smooth star formation until a single dom¬ 
inant burst at intermediate ages, Fornax has undergone at 
least two distinct bursts of significant star formation activ¬ 
ity. This result is in line with work by del Pino et al. (2013), 
whose analysis focused on stars relatively close to the core 
of Fornax, and measured two distinct periods of increased 
star formation, ~ 8 and ~ 10 Gyr ago. 

Using a matched filtering technique, we have mapped 
out the two red giant branch populations (called here B- 
RGB and R-RGB), as well as the horizontal branch (HB) 
and young main sequence (blue plume, BP). We find that the 
B-RGB - most likely corresponding to an older population - 
is significantly more spatially diffuse, extending all the way 
out to the Fornax tidal radius. Conversely, the R-RGB pop¬ 
ulation is more centrally focused, and displays a somewhat 
less elliptical shape than its B-RGB counterpart. The blue 
horizontal branch, which indicates an ancient stellar popula¬ 
tion, follows the same basic profile as the B-RGB. These two 
populations likely evolved together. Lastly, the blue plume 
stars are the most centrally located, and their density con¬ 
tours are significantly more circular than the other Fornax 
populations. 

We have also identified an additional overdensity of 
young stars, ~ 7' from the centre of Fornax. This overden¬ 
sity is very similar to the new feature detected in |de Bo^ 


et al. (2013), and was selected against the relatively dense 


central regions of Fornax by focussing on bright blue plume 
stars (above the horizontal branch). 

There is a need to pursue spectroscopic studies of For¬ 
nax in more detail, and to considerably larger radii. The 
intermediate age red giant population dominates in the cen¬ 
tral region, thus limiting our knowledge of the older popula¬ 
tion. In particular, by studying stars at these larger radii we 
can hope to get a firmer handle on the star formation his¬ 
tory of Fornax, and determine whether the counter-rotation 
- perhaps the strongest surviving argument for a significant 
merger in Fornax’s history - is real. 
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